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MODELO DO PÊNDULO SIMPLES 


A equação diferencial que descreve o movimento do pêndulo 
simples é dada pela segunda Lei de Newton, considerando a 
equação do momento desenvolvido. Pode ser dada por: 


mL2ô + bÅ + тд1ѕеп = и 


Sendo: 

m = Massa do corpo 

b = Constante de amortecimento dinámico 

L = Comprimento da barra ou raio do movimento 
mg сове 


g = gravidade 


u = acáo de controle ou momento aplicado 


MODELO DO PÊNDULO SIMPLES 


Considerando: 


т=а 1-0460-1, о modelo do pêndulo pode ser 
reescrito como, com ação de controle nula, visando encontrar 
o ponto de equilibrio: 


0 + 0 + ѕеп = 0. 


Dessa maneira, е preciso дезсгемег um modelo ет езрасо де 
estado para demais análises де estabilidade до ponto de 
equilibrio: 


та 
те созе 


(A) REPRESENTAÇÃO EM ESPAÇO DE 


ESTADOS 


А Ё | 28 я т х= х = Ө Ех = 0 
Dada a equação não-linear Ө + 0 + зепб = 0. Faça а substituição: 4 . =) ; | ж 
б == 9 == Ө Хо = 
Resolvendo a equação para x, utilizando as substituições acima, temos: Х2= — SİN Ху — X2 
P formal d d Е | х2 | 
ortanto, na forma de espaço de estados x = P m — sin X4 — x 


(B) DETERMINACÁO DOS PONTOS DE 


EOUILÍBRIO 


Nos pontos de eguilíbrio пао temos variacáo temporal das variáveis de estado ===) Х = [(фи)=0 


Assim, = [51] = 201 -0 


sin т = Э 


Solucáo: | ви : =) sin x =0 ==) в z 


sinx; — x, = 0 1= п 


сал = (0, 0) 


Portanto, temos dois pontos de equilíbrio: 
Ред? иж = (л, 0) 


(C) ESTABILIDADE DOS PONTOS DE 


EQUILÍBRIO 


Sera feita a linearização pelo |“ Metodo de Lyapunov зе (2210) шинээ) x = Ax+ Bu 
дА (Хед Чеа) ô fı (Хед, Неа) дЛ (Хед Неа) 
А Ер д] (Хед, иеа) дж 9x2 Ој (Xeg'Ueg) фи 
As matrizes А е B зао dadas por: А = — = = е В = == = 
Р дх Of2(XepUeg) 97 (Хеда Чед) ди ô f2 (Хед, Ига) 
дж 0x2 дм. 
Р ist |= | С | | 
ага о sistema X = | l-sinx — x, teremos: 


= 22 51 шан [0] 


(С) ESTABILIDADE DOS PONTOS DE 


EQUILÍBRIO 


Para o primeiro ponto de equilíbrio: Рен + (0, 0) А = б | ев = Д 
E -1 = 0 


E NE Autovalores complexos conjugados 
Os autovalores de A зао Ло = —0,5 +] £ E P a 
2 Parte real negativa 


Peq1= (0,0) е assintoticamente estável 


Se todos os autovalores de A tem parte real negativa, diz-se que Хед é assintoticamente estável, pois x > Хед 
quando t — %, e então o ponto de equilíbrio é chamado de sumidouro (sink) 


(C) ESTABILIDADE DOS PONTOS DE 


EQUILÍBRIO 


Para o segundo ponto de equilíbrio: Р, = (л, 0) E=) А = É 1 | е В = K 


A = 1,618 ==) шигээ reais 


Os autovalores de A sáo ћ = – 0,618 A > 0 


Peg2= (т, 0) е instável 


Um ponto de equilibrio é instável se um ou mais autovalores de A tiver parte real positiva. 


’ Se alguns dos autovalores, mas não todos, têm parte real positiva, enquanto o resto tem parte real negativa, Хед ё 
chamado de ponto de sela 


(D) RETRATO DE FASE PRÓXIMO AOS 
PONTOS DE EOUILÍBRIO 


LÁ 


OXIMO AOS 
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Posição angular (em graus) 


(F) SOLUÇÃO NO TEMPO (0;0) 


Posição angular da esfera (teta) 
100 
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Tempo (s) 


PĚNDULO: 2? MÉTODO DE LYAPUNOV 


” Define-se uma função de Lyapunov que descreva а 
energia do sistema. 


As energias envolvidas são a potencial e cinética, 
2 


5 o mv 
descritas respectivamente рог: mgh e ——. 
h = L(1 — со5(0)) = L(1 — cos(x1)) 


у = 0 = х, 


PĚNDULO: 2? MÉTODO DE LYAPUNOV 


Fazendo as devidas substituicoes, temos a eguacáo de Lyapunov: 
х2 
У(х) = (1 — cos(x)) + = 


Рага ит conjunto D = {x Є В:0 < ху < п} V(x) > 0, (х; = п não atende uma 
das definições de Lyapunov, pois У (х) = 0 para x = 0). 


О valor máximo de |cos(x,)| é |, portanto о valor mínimo para (1 — cos(x1)) é 0. 
Assim, os valores de x, não afetam o sinal da função V (x). 


PĚNDULO: 2? MÉTODO DE LYAPUNOV 


Мо ponto de equilíbrio Xeq = (0,0), а função V(x) = 0 


Calculando — se a derivada: 


ау (х) 
dt 


у= = sen(x1)xX1 + хх, = зеп(ж)хо + X (—x» — sen(x1)) =< 


0 


Portanto x = 0 е localmente estável no sentido de Lyapunov, mas sem garantia sobre 
estabilidade assintótica. 


PĚNDULO: 2? MÉTODO DE LYAPUNOV 


Uma das formas de testar se o sistema é assintoticamente estável е pelo 
principio de invariáncia de La Salle, sem exigir que a derivada da função de 
Lyapunov seja definida negativa. 


Seja l Е R, tal que а; СР, e um conjunto Е = (x E w,:V(x) = 0}. 


Ou seja, У (x) = —х2 = 0. Logo, Ух; Е D,V(x) = 0 зе, e somente se, x, = 0 


PĚNDULO: 2? MÉTODO DE LYAPUNOV 


A derivada da função V 
é semidefinida em ш. 
Logo toda solução do 
sistema iniciando em 
со ё atraída para о 
maior conjunto 
invariante em E. 


Como pode ser 
observado na imagem 
ao lado: 
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PĚNDULO: 2? MÉTODO DE LYAPUNOV 


Como a origem е o único conjunto 
invariante em E,a origem é um 
ponto de equilibrio 
assintoticamente estável localmente. 


De outro modo, uma vez que a 
derivada da função de Lyapunov é 
uma reta passando pelo ponto de 
equilibrio, a origem е o único 
conjunto invariante. 


